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Pluto’s surface is dominated by the huge, pear-shaped basin Sputnik Planitia. 
It appears to be of impact origin, but modelling has not yet explained its 
peculiar geometry. We propose an impact mechanism that reproduces 

its topographic shape while also explaining its alignment near the Pluto- 
Charon axis. Using three-dimensional hydrodynamic simulations to model 
realistic collisions, we provide a hypothesis that does not rely upon a cold, 
stiff crust atop a contrarily liquid ocean where a differentiated -730 km 
ice-rock impactor collides at low-velocity into a subsolidus Pluto-like 
target. The result is anew geologic region dominated by impactor material, 
namely a basin that (in a 30° collision) closely reproduces the morphology 
of Sputnik Planitia, and a captured rocky impactor core that has penetrated 
the ice to accrete as a substantial, strength-supported mascon. This provides 
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an alternative explanation for Sputnik Planitia’s equatorial alignment and 
illustrates a regime in which strength effects, in low-velocity collisions 
between trans-Neptunian objects, lead to impactor-dominated regions on 
the surface and at depth. 


In 2015, the New Horizons space probe revealed Pluto’s surface to be 
geologically complex’. It is dominated by an -1,200 km x 2,000 km, 
nitrogen-ice-filled basin named Sputnik Planitia (SP)' whose central 
plains are 3-4 km below its surroundings’. These values correspond 
to the surface of the nitrogen ice, whose thickness is unknown. The 
quasi-elliptical shape and mountainous rim resemble a degraded 
impact basin’, which has motivated several studies into its forma- 
tion using computational impact models*». These suggest that a 
400-km-diameter impactor could produce an ~800 km basin compa- 
rable in diameter to SP; however, they are two-dimensional simulations 
of head-on collisions and cannot reproduce the feature’s elongated 
morphology. 

Here we use three-dimensional impact simulations to consider 
more realistic impact angles and impactor structures and the com- 
plex material and dynamical interactions that ensue. We use the 
smoothed-particle hydrocode SPHLATCH (ref. 6), whichis the basis for 


awide range of planetary collisional studies’”°, and apply equations of 
state (EOSs) for rock (dunite) and water ice to represent the bodies. 
Strength has been shown to be important, even up to Mars-size tar- 
gets'°” so we include shear strength and plasticity in all our simulations 
using a Drucker-Prager-like yield criterion. Yield strength is a function 
of pressure and temperature, so is pronounced in cold materials and, 
thus, applicable to Pluto. The specific strength parameters follow 
those of previous impact studies” and are given in Extended Data 
Table 1. Further geophysical parameters, such as cohesion, tensile 
fracture and porosity, are important in collisions at the scale of aster- 
oids and comets” but are negligible in giant basin formation and are 
not included here. 

Pluto is represented as a differentiated body composed of a rocky 
core and a water-ice mantle. The former constitutes two thirds of 
the total mass in agreement with observations by New Horizons”. 
The physical state of the H,O mantle is uncertain, with recent studies 
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advocating for subsurface water that lies below a thick ice shell’*”’. One 
particular motivation for this scenario is the location of SP. The forma- 
tion ofa large topographic depression is expected to induce true polar 
wander, reorienting Pluto to position the depression closer to the near- 
est pole. Instead, SP is near Pluto’s equator and opposite its large, tidally 
locked, satellite Charon, implying the presence of amass concentration 
(‘mascon’)”’ rather than a deficit. The prevailing hypothesis is that a sub- 
surface ocean can explain such a mascon, as the excavation of SP caused 
an uplift of denser salty water from below*”’. The crater structure and 
subsurface ocean must persist into the present, because if the ocean 
had solidified or the ice crust had relaxed viscously, then the mascon 
would have vanished. Preventing the ocean from freezing necessitates 
factors suchas an exceptionally high ammonia concentration to lower 
the melting temperature of ice or an ever-present layer of gas hydrates 
that insulates the warm liquid ocean against the gelid ice shell”. 

To bracket the possible pre-impact conditions, our suite of simula- 
tions includes Pluto targets with solid rock cores (dunite EOS) and solid 
ice mantles between 70 and 250 K and cases with subsurface oceans 
with temperature profiles consistent with a water EOS. Impactors 
are modelled using the same materials, but their parameters are less 
constrained. Kuiper belt objects exhibit a wide range of diameters” ”, 
and despite observations of over 2,000 of these objects”, meaningful 
mass estimates are available for only ahandful”, making their densities 
and compositions quite uncertain. For this work, we consider solid 
impactors between 400 and 1,100 kmin diameter, with compositions 
ranging from rock toice and to in-between (5-66 wt% rock surrounded 
by ice). In the latter, multi-material cases, we model the impactor as 
fully differentiated, as is the case for Pluto. In reality, however, the 
degree of differentiation decreases with size”. Partial differentiation 
has even been suggested for bodies as large as Charon”. The impli- 
cations of this discrepancy are considered in the ‘Discussion’. In any 
case, the smaller size of the impactors implies that there is a cooler 
temperature than that of Pluto, and they are, therefore, modelled with 
ice-mantle temperatures of ~70 K. The precise interior profiles used 
forthe nominal impactor and target are given in Extended Data Figs.1 
and 2, respectively. Other target temperature profiles explored are 
given in Extended Data Fig. 3. 

We considered impacts with angles ranging from 0° to 45° 
and velocities of 1.0u,, to 1.4u,,,. Here, Vese represents the mutual 
escape speed for the two bodies (-1.2 kms‘), and is given by 


2G(Mimp+ Bags 
Vee =i] ee , where G denotes the gravitational constant, 
Timp ttar 


mis mass and ris radius, and the subscripts imp and tar signify the 
impactor and target, respectively. The impacts are slower than the 
speed of sound in geologic ice, -2-4 kms”, so shocks play a minor role. 
This velocity range is like that of Pluto—Charon formation scenarios”. 
If SP formed when Pluto was much closer to the Sun”, random encoun- 
ters would probably have been faster; in this case our scenarios might 
apply to a collision with a massive satellite of Pluto or to the last major 
event in Pluto’s original accretion, because the mass ratios, angles and 
velocities lie within the ‘merging’ impact regime**. 

The main impact sequence (shock wave traversal and dissipation 
and transient crater formation and collapse) occurs quickly, within the 
first -2 h after impact. The vast majority of displaced material settles 
soon thereafter, with only a small amount of ejecta on target-crossing 
orbits that will re-impact at later times. Our simulations were extended 
for 6 hafter impact, such that deformation has ceased. Any remaining 
lofted ejecta composes less than 0.5% of the total mass. 


Results 

Upon examination of the impact outcomes, distinct trends become 
evident. Figure 1 shows cross sections of the target for a variety of 
cases, from head-on at impact velocity 1.0v,,, (left) to oblique at impact 
velocity 1.4v,,, (right), for impactors of ice, rock and a differentiated 
ice-rock composition (top to bottom). In each case, a transient crater 


forms in Pluto’s mantle (larger for the larger-diameter ice impactor), 
which by t= 6 h after initial contact has collapsed. The ice impactor 
(Fig. la—c) experiences profound deformation upon colliding with 
Pluto’s ice mantle, spreading out to become a‘splat’ akin to impactor 
distributions proposed to explain the Moon’s striking nearside-far- 
side dichotomy”. The impacting material starts off less compressed 
and thus somewhat less dense, so remains mostly atop Pluto’s 
ice throughout the crater collapse. The rock impactor (Fig. 1d-f) is 
denser, stronger and more difficult to melt and is much more resist- 
ant to deformation and breaking up, so instead, it pierces Pluto’s 
ice mantle, remaining largely intact until it impacts Pluto’s rocky 
core. Then it spreads out into a hidden, deeply buried splat, a rocky 
mass atop the core-mantle boundary. 

Increasing the impact angle and velocity contribute a downrange 
motion, so that the distribution of the impactor material becomes 
asymmetric. In these relatively low-velocity (subsonic) collisions, most 
of the material of the ice impactor remains above the ice mantle of 
Pluto regardless of impact angle. For a larger impact angle and velocity, 
the break-up of the impactor and its continuing momentum can lead to 
a fraction of it going into orbit and forming a depositional ring along 
the impact plane where it re-impacts. 

Rock impactors also continue downrange in oblique collisions, 
but being denser, they plough through the icy mantle in a complex 
interaction. The details of these interactions depend on the material 
parameters of the colliding bodies and on the physics of mixing of 
dissimilar materials, with the added complexity that material remains 
cold and rigid well beyond the transient cavity. Generally speaking, for 
impact angles $30°, a rocky impactor comes toa halt downrange, and 
if it overcomes the strength of cold ice, it sinks towards the core (for 
example, Fig. 1e). For oblique collisions (230°) and for higher velocities, 
arocky impactor can slide and bounce downrange until it comes back 
out of the icy mantle before breaking up into a ring of impactor mate- 
rial that extends downrange from the impact site (for example, Fig. If). 

A differentiated impactor (Fig. 1g-i) results in scenarios that 
effectively combine the features of the ice and rock end members. As 
the ice is the first material to make contact with Pluto, it shields the 
impactor core from substantial deformation and heating, soit remains 
more intact than an impactor without an ice mantle. In the head-on 
case (Fig. 1g), the rocky core ends up at the core-mantle boundary, 
whereas the mantle ice fills the transient cavity. A notable distinction 
from the pure-ice scenario is the considerable amount of impactor 
ice pulled down to the core-mantle boundary by the impactor core. 
Consequently, the impact site comprises predominantly impactor 
material extending down to Pluto’s core. 

We find that these intermediate parameters, specifically with a 
core mass fraction of 5-30%, impact angle between 15° and 30° and 
an impactor diameter of around 700 km (comparable to the mass of 
asteroid Vesta), yield the most promising conditions for generating 
SP. An example in this range is shown in Fig. 2, which is our case that 
best fits the observations. The impactor first excavates the mantle of 
the impact site (Fig. 2a) to produce an elliptical transient crater. The 
core, shielded by the ice, stays largely intact as it is deflected through 
Pluto’s ice shell, ultimately halting while remaining embedded within 
the mantle (Fig. 2b). The transient crater becomes filled with infalling 
impactor ice, and the impactor core settles towards Pluto’s core to 
produce a mascon near the core-mantle boundary (Fig. 2c). 

The final outcome of these cases is a region of Pluto entirely domi- 
nated by impactor material with a shape in longitude-latitude space 
remarkably like that of SP, as shown in Fig. 3 (and Extended Data Fig. 4), 
which displays the final mass distribution of impactor material anda 
comparison with SP. The more circular shape in the north corresponds 
tothe point of collision, where the transient cavity opens and collapses. 
The pointed configuration in the southern region corresponds to the 
notable downrange movement of the impactor core through the target 
mantle, until it comes to rest atop the target core. 
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Fig. 1| Results of various impact simulations based on smoothed-particle 
hydrodynamics that highlight important regimes in the explored parameter 
space. a-i, Each plot shows the resulting body t= 6 h after initial contact asa 
cross-section of the impact plane with a thickness of 300 km. Colour denotes 

the composition and source parent body of the material, as labelled. The 
impactor masses are the same in each case, so that the ice impactor is the 

largest, approximately 730 km in diameter. The bodies are initially solid, and 


Target 


their temperature profiles correspond to the nominal case. The differentiated 
impactors havea core mass fraction of 15%. For an illustration of the initial 

set-up of these simulations, see Supplementary Fig. 1. a, Ice impactor, 1.0v,,., 0°. 
b, Ice impactor, 1.2v,,,, 30°. c, Ice impactor, 1.4V,,.,45°.d, Rock impactor, 1.0Vese 0°. 
e, Rock impactor, 1.2v,,., 30°. f, Rock impactor, 1.4v,.,,45°. g, Differentiated 
impactor, 1.0v,.,,, 0°. h, Differentiated impactor, 1.2v,,., 30°. i, Differentiated 
impactor, 1.4Vesc, 45°. 


The precise distribution of the impactor-dominated region is 
controlled by the impactor’s initial parameters. Its size determines the 
broad scale of the distribution, whereas its core mass fraction controls 
the breadth of the ‘tail’ of our teardrop-like distribution, owing to its 
dependence on core radius. The impact velocity also affects the scale 
of the transient crater, thus influencing the size of the elliptical region 
of impactor ice that forms at the immediate point of impact. However, 
its primary influence is on the length of the downrange tail. Higher 
velocities cause the core to travel a greater distance through the mantle 
before reaching a standstill and falling back towards the core-mantle 
boundary to create downrange disruption and surface deposition. 

Considering all these factors, we find that a 15% rock, ~730 km 
impactor striking Pluto at 30° with an impact velocity, Veo Of 1.2Vec¢ 
(-6 km s”) produces an impactor-dominated region at the impact 
site that best matches the planimetric morphology of SP (as in Fig. 3). 

In each of these cases, the rocky core of the impactor settles 
beneath the southernmost, narrow tail region of the impact site and 


becomes buried beneath the impactor mantle above. We propose that 
this material serves as the mascon responsible for SP’s current location 
on Pluto. The substantial mass excess provided by the dense rock is 
expected to drive SP’s position towards Pluto’s equator through true 
polar wander, ina similar fashion as described in previous studies”. 
In this scenario, we expect that the tail of SP (that is the region with 
the strongest mass excess) would align close to Pluto’s equator, as is 
observed today. 

Finally, the SP-shaped, impactor-dominated region must produce 
the basin associated with SP today. We propose that this could be due 
to a contrast in composition between the impactor and Pluto. If the 
impactor mantle was of a greater density than the Pluto mantle, the 
impactor-dominated region would apply a greater load on Pluto’s 
silicate core. This would cause it to sink to lower elevations than its sur- 
roundings as it moves towards isostasy, thus leading to the observed 
depression. As the load on Pluto’s core would be greatest in the region 
containing the impactor core, this region would initially be the deepest 
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Fig. 2 | Time series of a simulation based on smoothed-particle hydro- 
dynamics for an ~730-km-diameter impactor, 15% core mass fraction, 

Von = 1.20,,,and an impact angle of 30°. a-c, Slices through the impact plane. 
Colours indicate material composition and source parent body, with purple and 
green indicating impactor rock and ice, and blue and yellow indicating target rock 
andice, respectively, as in Fig. 1. The faded body shows the size and position of the 


1,000 km 


impactor at the moment of impact, with the arrow indicating the impact velocity. 
a, Shortly after impact, the transient crater is still present. The impactor core has 
impacted Pluto’s core and continues downrange. b, The transient SP crater has 
collapsed and been infilled with ice from the impactor. Beneath, the impactor 
core has nearly come to rest along the ice-rock boundary. c, State after impact, 
t=6h. The rocky impactor core is at rest under the narrow southern end of SP. 


depression within SP; however, N, ice would would rapidly accumulate 
in the basin”, smoothing out the surface and concealing the mantle’s 
previous topography. A diagram depicting the proposed internal 
structure of SP can be seen in Fig. 4. 

A contrast in density between the impactor and target material 
might be expected due to the different degrees of differentiation 
inside each body based on their very different masses. The impactor 
may, thus, have had a less-pure, more dusty ice mantle than the purer, 
and thus, lower-density, ice mantle of the pre-impact Pluto. Likewise, 
aless completely differentiated impactor core might have retained 
more lower-density constituents and been less dense, with a buoy- 
ancy against the higher-density rocky core of Pluto, further aiding its 
long-term stability at the core-mantle boundary. 


Discussion 

These results may explain SP’s unusual morphology and the proposed 
positive gravity contribution without the need for a subsurface ocean. 
The existence of an ocean would depend on the temperature profile 
of Pluto, which is not well constrained. Our nominal target has a solid 
ice mantle that starts at -130 K anda solid rock core at -195 K (Extended 
Data Fig. 2). This is consistent with ice-shell temperatures in previous 
SP formation models** and with temperature predictions from 
studies on the thermal evolution of Pluto, at least for early times 
(within roughly 250-500 Myr post-formation)**”” when impacts of this 
magnitude were probable”. The rock-core temperature is consistent 
with that from some previous models*”", but it has been suggested 
that after ~5OO Myr it may reach -500-700 K due to radiogenic heat- 
ing’’*'**, We, therefore, repeated our best-fitting impact simula- 
tion using a core temperature of ~500 K but found only negligible 
differences. 

To explain SP, the basin shape and profile as predicted by these 
models must be supported over geologic time. In cold solid materials, 
the basin structure (SP itself and the dense buried mascon provided 
by the impactor’s core) would be supported by strength. Isostatic 
readjustments would occur around these emplaced materials, with 
the mascon, thus, underlying the deepest region of the basin, and the 
impactor mantle (if it is slightly denser than Pluto’s mantle) would 
underlie a broader depression. N, ice would quickly accumulate inside 
the basin*”*’, potentially contributing further to the positive mass 
anomaly that we propose drove SP into its current position through 
true polar wander. The southernmost, narrow section of SP, directly 
above the mascon in our simulations, would be forced closest to the 
equator, matching observations. 

Temperature can affect our proposed scenario, even for solid 
targets, because the formation of SP depends on material strength 
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Fig. 3 | The distribution ofimpactor material after impact compared to 
Pluto’s observed elevation distribution. Top, distribution ofimpactor material 
after the nominal collision (same as Fig. 2), displayed as an equirectangular 
projection of thickness to a maximum depth of -150 km. The spherical reference 
frame has been chosen to best match the location and orientation of SP today. 
North of the equator, the distribution is nearly elliptical, corresponding to the 
impactor’s first point of contact and the subsequent transient crater. South of 
the equator lies the distribution’s pointed ‘tail’, which is caused by the impactor 
core sliding through the Plutonian mantle. Bottom, Pluto’s observed elevation 
distribution. SP is the feature in the centre. For a further comparison, see 
Extended Data Fig. 4. Figure reproduced with permission from ref. 2, Elsevier. 


(Extended Data Fig. 5), whichis weaker for higher temperature bodies 
and zero in an oceanic layer. We, therefore, tested the sensitivity of 
ourimpact mechanism to these effects by repeating the nominal impact 
with a broad range of target ice-mantle temperatures. For our cold- 
est targets, which start at ~70 K, the final material distribution of SP 
(Extended Data Fig. 6, top) looks very like our nominal case (Fig. 3). 
For somewhat warmer ice, ~190 K (Extended Data Fig. 6, middle), the 
final distribution also looks similar; however, a stronger uplift of the 
target mantle occurs beneath the impactor core. For even warmer but 
still solid ice, ~250 K, the impact induces substantial melting, resulting 
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Fig. 4 | Schematic diagram depicting the proposed geologic structure 
beneath SP. Colour signifies the composition, as indicated by the labels. The 
surface of SP corresponds with the surface of the N, ice. The ice mantle directly 
below SP is of impactor origin, according to the simulations. If the impactor 
mantle has asomewhat greater density than the ice mantle of Pluto, then the 
entire SP region will sink towards isostasy to create a regional depression. The 
impactor’s rock core ends up buried beneath the southernmost region of SP, 
where it remains. The higher density of this mascon, and the resulting deeper 
topography and reservoir of N, ice above it, leads to a positive gravity anomaly 
that establishes and maintains SP’s near-equatorial position through true 
polar wander. 


in much more mixing of Pluto’s original ice and water with that of the 
impactor. Consequently, the final distribution of impactor material 
becomes less distinct, having a more elliptical shape (Extended Data 
Fig. 6, bottom). These results indicate that our proposed impact mecha- 
nism is sensitive only to temperatures near (<50 K below) or above the 
solidus, at which the material strength is not adequately maintained. 
At lower temperatures, the impact scenario remains very much like 
our reference case. Therefore, differences between our simplified 
temperature profiles (Extended Data Fig. 2) and that of the real Pluto 
would not affect our results as long as Pluto’s interior was sufficiently 
subsolidus at the time of impact. 

We also considered targets that have a subsurface water ocean 
50, 100 or 150 km thick, as per earlier investigations’. The ocean 
thicknesses correspond toa temperature-pressure profile consistent 
with the water EOS, with a transition to a solid ice shell in the exterior 
(see Supplementary Fig. 2 for the exact interior profiles used). 

For a starting ocean 50 km thick, we observe that the impactor 
material distribution still resembles SP quite well, though with a shape 
that is more compact than for a solid target. The elliptical region 
produced by the impactor’s contact and compression into Pluto is 
particularly diminished (Fig. 5). This pattern persists for 100-km-thick 
and 150-km-thick oceans, which end up with an even smaller con- 
trast between the elliptical area around the region of contact and the 
tail sculpted by the continuing impactor core. These instances 
with thicker oceans bear minimal resemblance to SP’s morphology 
(Extended Data Fig. 7), and we, therefore, deem them unlikely. 

Overall, we find that a thin, target ocean (S50 km thick) is favoured 
among our candidate scenarios. This aligns with previous subsurface 
ocean thickness estimates derived from thermal evolution models of 
Pluto that assumed an ammonia content like that of other Kuiper belt 
objectss (<1%) and atypical reference viscosity for ice (10 Pa s)***"4*"*”, 
If the influence of ammonia proves negligible, Pluto might not possess 
a subsurface ocean at all*’*, in accordance with our nominal case. 

In our proposed scenario, the impact that produces SP leads to 
a band of re-impacting ejecta in the impact plane (Fig. 2). Reference? 
suggested that re-impacting SP material could be related to the origin 
of the complex ridge-trough system (RTS) that stretches north to 
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Fig. 5 | Equirectangular projections of the impactor material distribution 
after impact for each subsurface ocean thickness. Top, subsurface ocean 
thickness of 50 km. Middle, subsurface ocean thickness of 100 km. Bottom, 
subsurface ocean thickness of 150 km. These cases can also be viewed as an 
interior cross-section in Supplementary Fig. 3. 


south across Pluto’s surface and is interrupted at SP. SP and RTS have 
an -10° offset in longitude-latitude space, but this could possibly be 
accounted for by pre-impact rotation (see ref. 48 for asimilar scenario 
involving impact ejecta on Vesta). 

RTS has also been proposed to be a paleoequator’, meaning it 
could be linked to the true polar wander caused by SP (ref. 49). In this 
case, the impact plane of SP would be very close to the paleoequatorial 
plane, suggesting that the SP-forming impact may have contributed 
considerably to Pluto’s rotation. Clearly, future work including both 
pre-impact rotation and subsequent true polar wander calculations 
are necessary to further understand SP’s origin and its implications 
for other features such as RTS. 

Other follow-up work should include modelling the gravity field 
associated with the positive anomaly of the buried impactor core and 
exploring other possible compositions for the impactor and target. 
These results could then informa sophisticated geophysical model to 
investigate the long-term evolution of SP over geological timescales. 
This would be particularly important for understanding the proposed 
isostatic readjustment. 

Further impact studies with a wider range of impacts anda more 
systematic approach should also be performed to better constrain 
the parameter space of this regime. Higher velocity and more oblique 
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impacts will cause the impactor core to burst through the target 
planet’s surface and re-impact further downrange than the reference 
impactor (as in Fig. 1i); however, the limiting speed and angle for this 
scenario will change for different impactor sizes and compositions. 

Finally, we highlight that the impact mechanism outlined in this 
study is not limited to the formation of SP. Impacts between mostly 
solid bodies with masses up to that of Pluto and occurring at close tothe 
escape velocity have been shown to generate heterogeneous interior 
structures that may persist tothe current era. The proper interpretation 
of those structures may, thus, be crucial in discerning the properties 
(both interior and exterior) of the icy worlds that populate the cold 
outer fringes of planetary systems. 


Methods 
Smoothed-particle hydrodynamics 
The impact simulations were conducted with the three-dimensional 
smoothed-particle hydrocode SPHLATCH (ref. 6). This method repre- 
sents a continuum by a set of smoothed particles, each representing 
a fixed mass distributed across a sphere with radius 2h, where h is 
the smoothing length. The physical properties (for example density 
and internal energy) of an individual particle were computed using 
a kernel-based interpolation method that considers all neighbour- 
ing particles within the 2h-radius sphere. To ensure that there was a 
roughly constant neighbour count inside the smoothing sphere of 
each particle, h was variable, with regions of lower density having larger 
smoothing lengths. Self-gravity was included using a Barnes and Hut 
tree algorithm”. Distant particle contributions were calculated with 
the multipole approximation (up to the quadrupole moments). For an 
in-depth description of SPHLATCH, we direct interested readers to ref. 51. 
A total of 200,000 particles were employed in each simulation. 
To test for convergence, the best case was performed again using 
1,000,000 particles. This high-resolution run was utilized for the 
best case figures (Figs. 2 and 3); however, only negligible differences 
were observed between the two resolutions, and thus, we deem all our 
simulations adequately resolved. Each collision was simulated for 6 h 
after impact. 


Material properties 

Using the software ANEOS (refs. 52,53), an EOS was used to calcu- 
late pressure and temperature, with input parameters for dunite” 
and ice”***°. The latent heat of ice was accounted for within the ANEOS 
(refs. 57,58). The latent heat was not used for dunite, but this did not 
affect our results as all rocky material remained well below melting 
temperatures throughout the impact simulations. The material 
strength was represented by a Drucker-Prager friction yield crite- 
rion, which varies as a function of pressure and temperature’. The 
strength parameters are consistent with those from previous impact 
modelling studies”, as inferred by experiments conducted within 
the temperature and pressure range relevant to this work. A compre- 
hensive description of the strength model can be found in ref. 11 using 
the corrections of ref. 10. 

Extended Data Table 1 lists the input parameters used in our 
strength model. Due to the large, planetary scale of the impacts con- 
sidered in this study, we considered all material to be fully damaged. 
Dunite input values are from ref. 12. Our ice values are very close to 
those of ref. 15, which were determined by fitting to experimental 
results for quasi-static shear strength, friction and dynamic tensile 
strength. Other studies have used substantially higher yield parameters 
for ice”, Reference » adjusted the parameters from ref. 14 to better 
reflect the frigid temperatures of large icy satellites like Ganymede and 
Europa, which suit the conditions considered in this work. We make one 
modification, namely to the yield strength (Y) of damaged (fractured) 
ice, with the equation of ref. 62: 


Ya = min(Ya, Yan) (1) 


with 


Yar = HasP (2) 


and 


Yon =HahP + Youn» (3) 


where Y,,, and Y,, represent the yield strength of damaged ice at high 
and low pressures, respectively, P denotes pressure, Han and Ha, repre- 
sent the friction coefficients for damaged ice at high and low pressures, 
respectively, and Yo, represents the y intercept for an empirical fit to 
the experimental results at high pressure. 


Impactor distribution maps 

To calculate the regional distribution of impactor material onthe final 
Pluto target, we employed a spherical grid akin to that of ref. 10. An 
impactor mass fraction fimp was assigned for each particle. Impactor 
particles have fimp = 1, and the remaining particles that constitute the 
target have bimp = 0. This value can be calculated for any point in space 
with a smoothed-particle hydrodynamics sum: 


Bimp() a > Bimp W(X — Xp, hp), (4) 
b 


wherex denotes a given spatial point, b signifies a particular neighbour- 
ing particle and W represents the smoothing kernel, which is the 
cubic B-spline function. The spherical grid was then constructed with 
a resolution of 7 km in radius and — in latitude and longitude. imp 
was calculated for every grid point with equation (4). By summing 
all Bimp Values at each latitude and longitude coordinate down to a 
specified depth and multiplying this sum by the grid resolution in 
radius Ar, the total thickness of impactor material z,,,, was calculated 
at each coordinate: 


Zimp, jk = Ar >) Bimp,ijko (5) 


i=inin 


with subscripts i, j and k denoting the grid cell index for each radius, 
latitude and longitude, respectively, and imin denoting the index for the 
specified depth, which we take as ~150 km. 

To change the reference frame of our spherical coordinate system, 
we rotated the particle coordinates before constructing the spherical 
grid. This method was used on our best-fitting simulation case to match 
the location and orientation observed of SP, allowing for better visual 
comparison in Fig. 3. The reference frames of all subsequent impactor 
distribution maps (Fig. 5 and Extended Data Fig. 6) were also rotated 
by the same degree. 


Data availability 

The simulation data used to create all figures presented in this study has 
been published on Zenodo, https://doi.org/10.5281/zenodo0.10696642 
(ref. 63). 


Code availability 

A compiled version of the SPHLATCH code, along with the necessary 
input files, can be obtained from the corresponding author upon 
reasonable request. 
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Extended Data Fig. 1| Interior profiles for the nominal impactor. (15% core mass fraction). Colour signifies composition, where blue represents rock material and 
green represents ice. 


Nature Astronomy 


Article https://doi.org/10.1038/s41550-024-02248-1 


—————S eT 1.4 as 
= 3.0 ao oe 
a x 
S9 25 a 0.8 Ž 170 
x x o 
ma — 5 160 
> cD) 0.6 D 
£ 2.0 5 © 150 
e a 0.4 A 
(a9) v E 
A 1.5 à 0.2 £ 190 
130 
—s 0.0 NS 
0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 
Radius [x10 °? km] Radius [x10 ? km] Radius [x10 ? km] 


Extended Data Fig. 2 | Interior profiles for the nominal target body. Colour signifies composition, where blue represents rock material and green represents ice. 
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Extended Data Fig. 3 | Alternative temperature profiles investigated for energy (and thus temperature) properties of the ice-mantle are different to 
Pluto’s mantle. Alternative temperature profiles investigated (all but the the nominal case (Extended Data Figure 2); any differences in the rock-core 
nominal case) for Pluto’s mantle. Colour signifies composition, where blue properties, or in the density or pressure of the ice mantle, are negligible. 


represents rock material and green represents ice. Note that only the internal 
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Extended Data Fig. 4| The post-impact distribution of impactor material (Bottom) Pluto’s observed elevation distribution, where SP is the feature outlined 
compared to Pluto’s observed elevation distribution, with a dashed- by a blue dashed-line in the centre.The SP outline follows that of ref. 2; for a direct 
line outline highlighting the shape of SP. (Top) The distribution of comparison of these outlines, see Supplementary Figure 4. Figure reproduced 


impactor material after the nominal collision (same as Fig. 3), displayed as with permission from ref. 2, Elsevier. 
an equirectangular projection of thickness to a maximum depth of - 150 km. 
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Extended Data Fig. 5 | The nominal impact scenario as a time series without rock and ice, and blue and yellow indicating target rock and ice, respectively, as 
including strength. Slice through the impact plane. Colours indicate material in Fig. 1. Substantial deformation of Pluto’s core occurs throughout the impact 
composition and source parent body, with purple and green indicating impactor event, in stark contrast to our simulations that incorporate strength. 
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Extended Data Fig. 6 | The distribution of impactor material after the the top plot to the bottom (that is ~ 70 K (top), ~- 190 K (middle), - 250 K (bottom), 
nominal collision for all additional temperature profiles explored for a solid following the profiles in Extended Data Figure 3 from left to right). These cases 
Pluto mantle. The distribution is displayed as an equirectangular projection of can also be viewed as an interior cross-section in Supplementary Figure 5. 


thickness to a maximum depth of - 150 km. Mantle temperature increases from 
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Extended Data Fig. 7 | Equirectangular projection of the post-impact thickness (50 km, 100 km and 150 km from top to bottom), with a dashed-line 
impactor material distribution for each subsurface ocean thickness, with a outline highlighting the shape of SP. The SP outline follows that of ref. 2; for a 
dashed-line outline highlighting the shape of SP. Equirectangular projection direct comparison of these outlines, see Supplementary Figure 4. 


of the post-impact impactor material distribution for each subsurface ocean 
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Extended Data Table 1| Parameters used in the material strength model 


Parameter Description Value for ice Value for rock 
Yn Shear strength at P = o0 115 MPa 3.5 GPa 

Yo Cohesion 0 0 

Yo,h Y-intercept of high pressure model 8.3 MPa 0 

hi Coefficient of friction for intact material 2 2 

Ha, Coefficient of friction for damaged material at low pressure 0.55 0.8 

ld,h Coefficient of friction for damaged material at high pressure 0.2 0.8 

Tm Melting temperature 273.2 K 1500 K 

Ẹ Thermal softening parameter 1.84 1.2 

G Shear modulus 3.52 GPa 72 GPa 


The parameters for rock (dunite) are from ref. 12, and those for ice from ref. 15 and ref. 62. Some more recent studies”*”' have adopted substantially stronger ice models in the studied regime; 
however, little justification is provided for those values, and some appear to be erroneous (for example a damaged coefficient of friction, p4, of 0.8, which should be 0.2 according to their ref. 62). 
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